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ABSTRACT 

We report on interferometric and radial-velocity observations of the double- 
lined 51-d period binary (A) component of the quadruple pre-main sequence 
(PMS) system V773 Tau. With these observations we have estimated prelim- 
inary visual and physical orbits of the V773 Tau A subsystem. Among other 
parameters, our orbit model includes an inclination of 66.0 ± 2.4 deg, and allows 
us to infer the component dynamical masses and system distance. In particu- 
lar we find component masses of 1.54 ± 0.14 and 1.332 ± 0.097 M Q for the Aa 
(primary) and Ab (secondary) components respectively. 

Our modeling of the subsystem component spectral energy distributions finds 
temperatures and luminosities consistent with previous studies, and coupled with 
the component mass estimates allows for comparison with PMS stellar models 
in the intermediate-mass range. We compare V773 Tau A component properties 
with several popular solar-composition models for intermediate-mass PMS stars. 
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All models predict masses consistent to within 2-sigma of the dynamically deter- 
mined values, though some models predict values that are more consistent than 
others. 

Subject headings: binaries: spectroscopic — stars: fundamental parameters - 
stars: pre-main sequence — stars: individual (V773 Tau) 

1. Introduction 

Accurate determinations of the physical properties of stars (e.g. mass, radius, temper- 
ature, luminosity, elemental abundances, rotation, etc.) provide fundamental tests of stellar 
structure and evolution models. Among the areas where our understanding of stellar struc- 
ture is most uncertain is in pre-main sequence (PMS) stars, particularly for low-mass systems 
(see Palla & Stahler 2001; Hillenbrand & White 2004; Mathieu et al 2007, and references 
for summaries). Previously we have reported results from our program of studying PMS 
binary systems with the intent of establishing dynamical mass and luminosity constraints 
(Boden et al. 2005b). Here we report on observations of the double-lined spectroscopic binary 
V773 Tau A, part of the PMS quadruple system V773 Tau. 

V773 Tau (HDE 283447, HBC 367) is among the most interesting and well-studied sys- 
tems in the Taurus-Auriga star forming region. V773 Tau was first identified as a T Tauri 
star by Rydgren et al (1976) based on Ha and Ca II H and K emission, high lithium abun- 
dance, photometric variability, and K2 spectral type. The object presented an enigmatic 
mixture of classical (CTTS) and weak-lined T Tauri (WTTS) characteristics until it became 
apparent there are multiple components; V773 Tau was resolved as a visual binary inde- 
pendently by Ghez et al (1993) and Leinert et al. (1993) (with components designated here 
as A and B). Martin et al (1994) suggested, and Welty (1995, herein W1995) established 
A as a short-period (51 d) double-lined spectroscopic binary. In high-angular resolution 
studies Duchene et al (2003, herein D2003) and Woitas (2003) independently identified an 
additional, "infrared" component in the system (herein designated C - note D2003 use an 
alternate component notation), making V773 Tau an apparent compact quadruple system 
with four stars within roughly 100 AU 1 . Using spatially-resolved photometry, D2003 find 



1 There are at least three different component notations used in describing the V773 Tau system. Herein 
we have described the 51-d SB2 subsystem as "A", with stellar components Aa and Ab; this is consistent 
with nomenclature used in Massi et al (2006) and Woitas (2003). Conversely, D2003 refer to the 51-d SB2 
subsystem as "AB". Finally, based on arguments made by Hartigan et al (1994) on the "likely" association 
of V773 Tau and FM Tau, the Washington Double Star Catalog (Mason et al 2001) describes the four-star 
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the A, B, and C components to have distinctly different spectral energy distributions. D2003 
conclude that V773 Tau contains stars at very different PMS evolutionary states: the opti- 
cally brighter A stars are WTTS in character, while the circumstellar material inferred from 
the significant IR excess is associated with the B and C components. 

Further, the V773 Tau A subsystem is a bright, polarized, and highly- variable radio 
source. Kutner et al (1986) discovered radio continuum emission from the system, appearing 
as the most luminous Taurus-Auriga T Tauri star in the survey by O'Neal et al (1990). 
Feigelson et al (1994) published a comprehensive multiwavelength study of the V773 Tau 
system, observing both polarized and flaring radio emission. W1995 argued that both A 
components are chromospherically active based on variations in their optical spectrum and 
estimates of relatively rapid rotation (v sini ~ 40 km s" 1 ); this point was confirmed by 
Phillips et al (1996, herein P1996), who resolved the radio emission into a clear binary 
morphology apparently corresponding to the two A components. Massi et al (2002) identified 
an increase in flare activity around A periastron, and Massi et al (2006) studied one such 
flare in detail; more frequent and powerful radio flaring at periastron is apparently due to 
interaction of component magnetic fields while in close proximity. Finally, Lestrade et al 
(1999, herein LI 999) used Very Long Baseline Interferometry (VLBI) astrometry to measure 
the trigonometric parallax for V773 Tau A as 6.74 ± 0.25 mas, corresponding to a system 
distance of 148.4 ± 5.5 pc. 

Here we report on observations of the V773 Tau A subsystem made with the Keck 
Interferometer (KI, see Colavita et al 2003), the Very Long Baseline Array (VLBA, see 
Napier et al 1994; Zensus et al. 1995), and spectroscopic radial velocities from the Harvard- 
Smithsonian Center for Astrophysics (CfA) "Digital Speedometers" (Latham 1992). These 
observations resolve V773 Tau A and allow us to estimate the visual and physical (i.e. three- 
dimensional) orbits, and determine the component dynamical masses and system distance. 
We describe these orbit models as preliminary because they weakly depend on ongoing 
observations and modeling of the larger V773 Tau system. Combining derived dynamical 
masses and system distance with component radiometric modeling, we estimate component 
temperatures and luminosities, allowing us to make direct comparisons with a number of 
popular solar-composition PMS stellar models. 



complex as "A" (with FM Tau as "B"), and the 51-d SB2 as "Aa" with components "Aal" and "Aa2". 
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2. Observations and Orbital Solution 

KI Observations The KI interferometric observable used for these measurements is the 
fringe contrast or visibility (squared, V 2 ) of an observed brightness distribution on the sky. 
Orbital analysis methods for such wide-band V 2 observations are discussed in Boden et al. 
(2000) and not repeated here. V773 Tau A was observed by KI in K-band (A ~ 2.2/zm) 
on five nights between 2002 Oct 24 and 2006 Dec 8, a data set spanning roughly four years 
and 30 orbital periods. V773 Tau A and calibration objects were typically observed multiple 
times during each of these nights, and each observation (scan) was approximately 130 sec 
long. For each interferometric scan we computed a mean (incoherent wide-band; Colavita 
1999; Colavita et al 2003) V 2 value from the scan data, and the error in the V 2 estimate is 
inferred from the rms internal scatter. For this analysis we have used HD 27741 (GO V) and 
HD 28483 (F6 V) as calibration objects. Calibrating our interferometric data with respect 
to these objects results in 27 calibrated visibility scans on V773 Tau A over the five epochs. 
Because the AO-corrected Keck Telescopes separately resolve the V773 Tau A-B-C complex 
(D2003), and the KI beam combiner is fed by single-mode fiber, no significant light from 
the B or C components enters the fringe camera while observing A (the if-band projected 
FWHM size of the KI single-mode fiber is 50 mas, so the B component is approximately 
five e-folding lengths away from A at the epoch of our observations; Colavita et al 2003, 
D2003). Consequently no special provisions for extra light are necessary in processing KI 
observations of V773 Tau A. 

Our KI V 2 observations are summarized in Table 1. The V 2 data set is depicted in 
Figure 1, along with V 2 predictions from our "Hybrid" orbit model (discussed below). 

VLBI Observations P1996 observed V773 Tau A with the VLBA, and with similar spa- 
tial frequencies as our KI observations (~ 50 MA) resolved the typical 2 mas separation 
between the A components. P1996 Figure 6 depicts a VLBA contour map of V773 Tau A 
at a particular epoch (1992 Sept 11), and clearly shows a binary source structure. P1996 
interpreted this binary morphology as radio flare emission from the two subsystem compo- 
nents. P1996 argued that the radio emission is coincident with the stars (to within a "few" 
stellar diameters), so a component separation interpreted as at-epoch relative astrometry 
was estimated by P1996. 

A team led by one of us (M.M.) has studied the physics of V773 Tau A radio flares. In 
addition to previous results (e.g. Massi et al 2002, 2006), V773 Tau A was observed using 
the VLBA in coordination with the Effelsberg 100 m radio telescope (herein VLBA+EB). 
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Fig. 1. — KI Data/Model Comparisons for V773 Tau A V 2 Observations. Here we give 
comparisons for our five epochs of KI V 2 observations and V 2 predictions from our "Hybrid" 
orbital model (Table 4). In each case the data and model are shown. 
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These 8.4-GHz (X-band) observations were conducted around apastron on seven consecutive 
days (2004 March 11-17 inclusive). As in the case of P1996, analysis of these observations 
revealed what appears to be both A components simultaneously exhibiting radio flares in 
most of these observations. However some of these "double-flare" observations show complex 
emission morphologies, complicating estimation of the binary relative separation. Reviewing 
the maps from these seven epochs, we decided to use astrometry from observations where 
there was significant (i.e. that the peak emission was greater than six sigma above the rms 
noise level) emission from both components, and the component structures were largely 
unresolved compared to the VLBA synthetic beam. These criteria resulted in our retaining 
three of the seven 2004 epochs for orbital analysis purposes; the astrometry resulting from 
these observations is summarized in Table 2. A sample image constructed from one of the 
three astrometric epochs (2004 March 15) is shown in Figure 2. The astrometric results of 
both the P1996 and these VLBI observations are depicted in Figure 3; these observations 
agree well with our orbit model for V773 Tau A described below. The radio morphologies 
seen in the full VLBI+EB data set will analyzed in a forthcoming publication (Massi et al 
2007, in prep). 

RV Observations W1995 established V773 Tau A as a double-lined spectroscopic bi- 
nary, and published an orbital solution. RV monitoring of V773 Tau with the CfA "Digital 
Speedometers" (Latham 1992) was carried out at the Oak Ridge Observatory (Harvard, 
Massachusetts), the F.L. Whipple Observatory (Mount Hopkins, Arizona), and the MMT 
Observatory (also on Mount Hopkins, Arizona). The nearly identical instruments at these 
facilities record a single echelle order and provide a wavelength coverage of 45 A centered at 
5187 A, with a resolving power of A/AA 35,000. Wavelength calibration was achieved in 
the usual manner through exposures of a Th-Ar lamp taken before and after each science 
exposure. The velocity zero point was monitored by means of exposures of the dusk and 
dawn sky, and small run-to-run corrections were applied to the velocities derived below as 
described by Stefanik et al. (1999). 

Our spectroscopic observations of V773 Tau began roughly at the same time as the 
W1995 data set and have continued to the present. In total we have collected 52 spectra 
spanning the period 1991 Sep 25 to 2007 Jan 7, with signal-to-noise ratios of 16-35 per 
resolution element of 8.5 kms -1 . The Radial Velocities (RVs) derived from these observations 
are summarized in Table 3. RVs were derived using the two-dimensional cross-correlation 
technique TODCOR (Zucker & Mazeh 1994), with templates adopted from a large library 
of calculated spectra based on Kurucz model atmospheres (see also Nordstrom et al. 1994; 
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Fig. 2. — Sample VLB A Map from 2004 March 15. Here we give a contour map of the total 
intensity continuum emission toward V773 Tau A from VLBA+EB observations of 2004 
March 15. Contours are at -3, 3, 4.5, 6, 7.5, 9, 12, and 15 a of the rms noise in the image 
(0.1 mJy/beam). The FWHM beam size (lower left) is 1.36x0.50 mas at a position angle 
of —11.6°. The peak flux density in the image is 0.55 mJy/beam. V773 Tau A is offset in 
these data and image as a calibrator (J0403+2600, not shown) served as the phase center. 
Heavy crosses are rendered at the predicted stellar separation based on our "V 2 & RV" orbit 
from Table 4 (i.e. no radio data included in the orbit modeling), registered to the center 
of the apparent primary position. The sizes of the crosses indicate the uncertainty in the 
separation estimate (~ 0.1 mas per axis). 
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Fig. 3. — Orbit of V773 Tau A as derived in our "Hybrid" solution (Table 4). Upper Panel: 
the relative visual orbit model of V773 Tau A is shown, with the primary and secondary 
objects rendered at T (periastron). The specific epochs where we have KI V 2 phase coverage 
are indicated on the relative orbit (they are not separation vector estimates), and the VLBI 
separations from P1996 (in blue) and this work (in red) are also shown. Component diameter 
values are estimated and rendered to scale. Lower Panel: the double-lined radial velocity 
orbit model and data described here. 
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Latham et al. 2002). The best-matching templates were chosen by running large grids 
of cross-correlations as described by Torres et al. (2003), which also yielded estimates of 
some of the stellar parameters. For an assumed Taurus solar metallicity (see Padgett 1996) 
and surface gravities appropriate for these stars (log g ~ 4.0; see Table 5), we derived 
effective temperatures of 4900 ± 150 K and 4740 ± 200 K for Aa and Ab, respectively. The 
adopted errors are conservative estimates to account for the abundance uncertainty. Formal 
errors are 50 K smaller, and are based on the scatter from the different spectra. Using 
the same methods we estimated the projected rotational velocity for the primary star to be 
vsini = 38 ± 4 kms -1 , in good agreement with the determination by W1995 (41.4 kms -1 ; 
no uncertainty reported). We were not able to determine the vsini for the fainter secondary, 
possibly due to the relatively low signal-to-noise ratios of our spectra. The templates for our 
radial-velocity determinations were selected to have temperatures of 5000 K and 4750 K (the 
nearest values in our template library) for Aa and Ab, and rotational velocities of 40 kms" 1 
for both stars. The light ratio between the secondary and the primary was determined with 
TODCOR by leaving it as an additional free parameter during the velocity determinations, 
as described by Zucker & Mazeh (1994). We obtained iAb/^Aa = 0.37 ± 0.03 at a mean 
wavelength of 5187 A. 

The presence of two additional spatially unresolved components in the V773 Tau system 
may in principle affect both our determination of the stellar parameters of Aa and Ab and 
their radial velocities. While the infrared companion V773 Tau C (D2003, Woitas 2003) is 
much too faint to cause any contamination, V773 Tau B is estimated to be ~ 1.7 magnitudes 
fainter than the spectroscopic binary in the optical (D2003), and could possibly have an 
effect. We attempted to detect the lines of this component in our spectra using an extension 
of TODCOR to three dimensions (see Zucker et al. 1995), but were unsuccessful. On this 
basis we conclude that it does not significantly affect our results above. However, its presence 
does complicate the analysis of the RV data for the spectroscopic binary in that the motion 
in the A-B orbit needs to be taken into account. We discuss this further below, and in a 
forthcoming publication (Torres et al 2007, in prep). 

Orbital Solution We have combined the astrometric (KI, P1996, VLBA+EB) and radial 
velocity data described above to estimate the visual and physical (i.e. three-dimensional) 
orbits for V773 Tau A. In particular, methods for modeling the wide-band visibility data 
used in this analysis are described in Boden et al. (2000), and not repeated here. Our 
"Hybrid" orbital solution (Table 4) is depicted in Figure 3. The upper panel depicts the 
relative visual orbit model, with the primary (Aa) component rendered at the origin, and 
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the secondary (Ab) component rendered at periastron. The KI V 2 data phase coverage 
is indicated on the visual orbit with points (they are not separation vectors); the phase 
coverage of the V 2 data is sparse relative to other similar analyses (e.g. Boden et al. 2000), 
but phase coverage provided by VLBI observations from P1996 and this work complements 
the V 2 coverage. Further, the incorporation of the VLBI separations also breaks the 180°- 
Q degeneracy inherent in V 2 -onlj analyses (e.g. Boden et al. 2000, 2005a), so the A-orbit 
rendered in Figure 3 is indeed as it appears on the sky. The apparent size of the V773 Tau A 
components are estimated (§ 3.1, Table 5) and rendered to scale; these same diameter values 
are used in the orbital modeling. The lower panel depicts radial velocity curves from the 
"Hybrid" orbit model and radial velocity data described above. 

V773 Tau A presents two significant challenges for RV orbit estimation. The first is 
that both components are relatively rapid rotators: W1995 found the v sin i for both stars 
to be around 40 km s -1 ; this rotational broadening is seen in our spectra as well, and it 
significantly degrades the precision of the velocities. Secondly, and most notably, is the A- 
barycenter motion that results from the A-B orbit. D2003 Figure 4 depicts the significant 
A-B orbital motion over the time since the detection of V773 Tau B, and our RV data set 
(starting in 1991) not only spans essentially all of this period, but includes the A-B periastron 
passage expected in 1996 (D2003). It is therefore important to account for the A-B orbital 
motion in the modeling of the A component RVs. Modeling of the A-B orbit would not have 
been possible at the time of W1995, but becomes feasible with the present data set. 

A separate RV-only orbital solution was carried out in which, in addition to solving 
for preliminary elements of the A orbit, we modeled the Keplerian motion of A around the 
A-B barycenter. In this analysis we held fixed some of the elements of the outer orbit to 
the values given by D2003 (specifically, the period, eccentricity, longitude of periastron, and 
time of periastron passage) and we solved for the velocity semi-amplitude of the barycenter 
of A [Ka — 6.11 ± 0.34 km s _1 ). This motion was then removed from our original velocities 
prior to combining them with the KI and VLBI data. A more complete treatment would 
incorporate all the astrometric measurements for the A-B orbit (e.g., D2003) along with the 
velocities and astrometry for the inner subsystem in a full solution solving for both orbits 
simultaneously. Such a study is beyond the scope of the present work, but will be described 
in a forthcoming publication (Torres et al 2007, in prep). 

There are significant differences between several of our orbital parameters for V773 Tau A 
and the corresponding values from W1995, most notably in the period and K^ a (see Table 4). 
We attribute these differences to a combination of at least three effects: i) limited phase cov- 
erage of the A orbit by W1995, particularly near A apastron; ii) possible biases in the W1995 
RVs due to line blending, exacerbated by the rapid rotation of Aa and Ab along with the 
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lower spectral resolution of that study compared to ours; and Hi) residual effects from the 
un-modeled motion in the outer orbit in W1995's solution. 

PI 996 argue that their resolved VLBI image of V773 Tau A measures the relative 
component positions to within a "few" stellar diameters (the primary apparent diameter is 
approximately 0.15 mas; see § 3.1). However, it is important to recognize that VLBI and KI 
V 2 data measure different emission mechanisms, and further, one would not expect individual 
radio flares to be centered on the stellar photospheres (measured by the KI V 2 data). This 
potential photosphere/radio offset motivated us to approach integrating the radio data into 
the orbit model cautiously. In our initial screening of the VLBI relative astrometry we 
compared the radio observations with position predictions from a V 2 and RV orbit model 
(Table 4). We found the model reasonably reproduced the observed radio separations with 
a typical residual of 0.2 - 0.3 mas 4 R*), slightly larger than the formal errors expected 
from the VLBI synthsized beam and the snr in the image, but in good agreement with the 
expected IR/radio offsets. The agreement we see between the radio separations and the V 2 
& RV orbit model makes it clear that the VLBI data do contain useful information on the 
component separations. We therefore felt confident in integrating the KI V 2 and four VLBI 
observations (including the PI 996 epoch) into our orbit model estimation. An analysis of 
the radio mophologies seen in the VLBA+EB observations will be presented in a future 
publication (Massi et al 2007, in prep). 

Figure 1 depicts direct comparisons between our KI V 2 observations and predictions 
from our V773 Tau A "Hybrid" orbit model from Table 4 (the five nights of data are each 
rendered in separate sub-panels). The model is seen to be in good agreement with the KI 
data. Further, Figure 3 shows a direct comparison between the VLBA separations and the 
"Hybrid" visual orbit model. Again, the body of the VLBI data are in good agreement with 
the orbit model, and the VLBI phase coverage complements the phase coverage provided by 
the KI V 2 data. 

Orbital parameter estimates for the V773 Tau A subsystem are summarized in Table 4. 
Included for relative comparison are the orbital parameters determined by W1995, our RV 
double-lined orbit ("RV"), the KI V 2 & RV solution ("V 2 & RV"), and the full integrated 
solution including the RV, V 2 , and VLBA data ("Hybrid"). The "Hybrid" designation 
for the composite orbit model recognizes its unique inclusion of both near-IR and radio 
interferometry data; this model is our favored orbit solution and is used in the remaining 
analysis. Again note that the A barycenter velocity is artificially near zero in our orbital 
solutions as the RV measurements have been corrected for our model of the A-B orbital 
motion. 
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Table 4. V773 Tau A Orbital Parameters. 



Orbital 


W1995 




This Work 




Parameter 




RV 


V 2 & RV 


Hybrid 


Period (d) 


51.075 ± 0.018 


51.1055 ± 0.0023 


51.1033 ± 0.0022 


51.1039 ± 0.0021 


T (MJD) 


49330.44 ± 0.56 


53059.87 ± 0.36 


53059.50 ± 0.35 


53059.73 ± 0.33 


e 


0.267 ± 0.016 


0.275 ± 0.012 


0.2708 ± 0.0094 


0.2717 ± 0.0085 


K Aa (km s _1 ) 


32.6 ± 0.7 


35.95 ± 0.53 


35.96 ± 0.53 


35.90 ± 0.53 


K Ab (km s -1 ) 


43.1 ± 2.1 


41.6 ± 1.4 


41.5 ± 1.4 


41.5 ± 1.4 


7 (km s -1 ) 


13.8 ± 0.9 


0.00 ± 0.31 


0.05 ± 0.31 


0.02 ± 0.32 


^Aa (deg) 


10.5 ± 5.1 


5.1 ± 2.6 


3.4 ± 2.5 


4.6 ± 2.4 


ft (deg) 






61.5 ± 2.5 


63.5 ± 1.7 


% (deg) 






65.9 ± 2.8 


66.0 ± 2.4 


a (mas) 






2.777 ± 0.056 


2.811 ± 0.047 


A K (mag) 






0.559 ± 0.082 


0.551 ± 0.078 



Note. — Summarized here are orbital parameters for the V773 Tau A subsystem as deter- 
mined by W1995 and present results. We give three separate fits to our data sets: RV only, 
Kl V 2 integrated with RV ("V 2 & RV"), and Kl V 2 , VLB1, and RV ("Hybrid") . Barycenter 
velocity (7) estimates (shown in italics) are artificially near zero because we have estimated 
and removed a model of the A-B orbital motion (see discussion in text). VL is quoted in a 
position angle convention. 
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3. V773 Tau A physical properties 

The orbital parameters from Table 4 allow us to directly compute many of the physical 
properties of the V773 Tau A subsystem and its components. Physical parameters derived 
from our "Hybrid" integrated visual/spectroscopic orbit are summarized in Table 5. The 
dynamical masses resulting from our orbit model for the two A components are 1.54 ± 0.14 
and 1.332 ± 0.097 M© for Aa and Ab respectively. This primary mass is in acceptable 
agreement with inferences by White & Ghez (2001, herein WG2001) and Palla & Stahler 
(2001) based on radiometric properties (i.e. luminosity and effective temperaure). Because 
it bears on the A-B orbit analysis (Torres et al 2007, in prep) we estimate the A subsystem 
mass as 2.87 ± 0.18 M . 

L1999 estimated a V773 Tau system distance of 148.4 ± 5.5 pc. The distance determi- 
nation to V773 Tau A derived from our visual and physical orbital solution is 136.2 ± 3.7 
pc, corresponding to an orbital parallax of 7.34 ± 0.20 mas, and consistent with the L1999 
result at 8.1% and 1.9-sigma. L1999 describe their distance estimate error as "formal", and 
it apparently did not account for the orbital motion of A in their data reduction. With these 
caveats we believe the L1999 distance estimate and ours are in acceptable agreement. 

At the distance of V773 Tau, neither of the A subsystem components are significantly 
resolved by the KI K-band fringe spacing, and we must resort to model diameters for the 
components. We have estimated the V773 Tau A component apparent diameters through 
spectral energy distribution modeling; details of the spectral energy distribution modeling 
are given in §3.1. We find apparent diameters of 0.152 ± 0.014 and 0.119 ± 0.011 mas for 
the Aa and Ab components respectively. With our system distance estimate these estimated 
diameters correspond to physical radii of 2.22 ± 0.20 and 1.74 ± 0.19 R Q , and (combined 
with the mass estimates) log surface gravities of 3.930 ± 0.094 and 4.081 ± 0.088 for the 
Aa and Ab components respectively. Similarly we have estimated absolute magnitudes in 
V and K, and a V — K color index for each of the components (Table 5), but these are of 
limited precision as they are dominated by the uncertainty in the extinction to the system 
(§3.1). 



3.1. Spectral Energy Distribution Modeling 

Because interferometric observations potentially resolve the stellar components in a 
binary system, we always construct spectral energy distribution (SED) models for binary 
systems to estimate a priori apparent diameters. Sample results of our SED modeling of 
V773 Tau A are shown in Figure 4. Our modeling is based on visible and near-IR photometry 
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Table 5. V773 Tau A Physical Parameters. 



J. ii y oxv^ciix 


Primary 


Secondary 


Parameter 


(Aa) 


(Ab) ' 


a (AU) 


0.1777 ± 0.0042 


0.2053 ± 0.0070 


Mass (M ) 


1.54 ± 0.14 


1.332 ± 0.097 


System Mass (M ) 


2.87 ± 0.18 


System Distance (pc) 


136.2 ± 3.7 


ir orb (mas) 


7.34 ± 0.20 


T cS (K) 


4900 ± 150 


4740 ± 200 


Model Diameter (mas) 


0.152 ± 0.014 


0.119 ± 0.011 


Bolometric Flux (10 _9 erg cm" 2 s" 1 ) 


4.44 ± 0.40 


2.38 ± 0.15 


Luminosity (Lq) 


2.56 ± 0.35 


1.37 ± 0.15 


Radius (R©) 


2.22 ± 0.20 


1.74 ± 0.19 


logs' 


3.930 ± 0.094 


4.081 ± 0.088 


M^ (mag) 


1.75 ± 0.12 


2.30 ± 0.13 


M v (mag) 


3.90 ± 0.16 


4.98 ± 0.16 


V-K (mag) 


2.15 ± 0.14 


2.68 ± 0.14 



Note. — Summarized here are physical (and related) parameters for the 
V773 Tau A subsystem as derived primarily from the "Hybrid" orbit solution 
(Table 4) and radiometric modeling (§ 3.1). 
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on A from WG2001 and D2003 respectively (summarized in Table 6), and component flux 
ratios from W1995 and our visible spectroscopy and near-IR KI V 2 data (discussed above). 
Using a custom two-component SED modeling code, we have modeled the A subsystem 
flux and ratios using solar-abundance, intermediate surface gravity SED templates from 
Kurucz (2001) and Lejeune et al (1997, 1998). We model the flux and ratio data with a 
large grid of SED templates for each component ranging in effective temperature (spectral 
type) from 4250 - 5500 K, surface gravity of log g = 4.0, and including extinction as a free 
variable in the parameter estimation. Constrained by the ratio data, this SED modeling 
preferred component temperature differences over a reasonably narrow range (A T c s ~ 
150-250 K) consistent with our spectroscopic study. But because the component colors are 
somewhat degenerate with reddening from dust extinction we found acceptable solutions over 
a relatively broad range (~ 500 K) of component temperatures. We have therefore adopted 
templates corresponding to the spectroscopic temperatures (4900 and 4740 K for Aa and 
Ab respectively). These temperatures acceptably agree with the composite temperature 
estimated by WG2001, but the hotter primary in our modeling prefers a slightly higher 
extinction value (A y = 1.80 ± 0.15 vs 1.39 ± 0.17 from WG2001). The resulting component 
bolometric fluxes and luminosities (the second based on our distance determination) are 
summarized in Table 5. In particular we find luminosities of 2.56 ± 0.35 and 1.37 ± 0.15 L Q 
for Aa and Ab respectively; these values are in good agreement with the range of previous 
component luminosity estimates from W1995, Ghez et al (1997), and D2003. 



4. Discussion 

4.1. Comparison with stellar evolution models 

Hillenbrand & White (2004) and Mathieu et al (2007) have recently summarized com- 
parisons between PMS stellar models and PMS stars with dynamical mass determinations. 
In particular Hillenbrand & White (2004) has shown that PMS models are in good agree- 
ment with dynamically-determined masses above 1.2 M . However, below 1.2 M Q the exist- 
ing models do a poorer job of matching observed component properties, tending to predict 
hotter and more luminous stars for a given mass. 

Figure 5 presents theoretical HR-diagram (luminosity/T c g) comparisons between our 
determinations of V773 Tau A component radiometric properties and solar-composition PMS 
models from Siess et al (2000), Yonsei-Yale (Demarque et al 2004, herein Y 2 ), D'Antona & 
Mazzitelli (1997), and Palla & Stahler (1999). (Popular models by Baraffe et al. 1998, were 
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Fig. 4. — Sample Spectral Energy Distribution Model for V773 Tau A. Here SED templates 
from Lejeune et al (1997, 1998) have been used to simultaneously model published flux 
measurements (from WG2001 & D2003) and flux ratio estimates (Table 4; bottom) from our 
KI V 2 and spectroscopic measurements. 
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excluded only because they do not cover the mass range for these stars.) These models span 
a wide range of physical ingredients and assumptions, some more realistic than others for the 
particular mass range of interest. In these figures evolutionary tracks at fixed mass between 
1.2 and 1.8 M (in steps of 0.1 M as available from the relevant models) are shown in solid 
lines, with heavy lines indicating tracks at masses close to dynamical masses from Table 5. 
Superimposed in dotted lines are isochrones at a relevant range of ages (1, 3, 5, 7, and 10 
Myrs as available from the relevant models). Finally, the Y 2 models support computing 
tracks at user-defined masses; tracks at the estimated component masses from Table 5 are 
shown for Aa (1.54 M , blue) and Ab (1.332 M , red). 

As expected from Hillenbrand & White (2004), all the surveyed PMS models are in rea- 
sonable (e.g. two-sigma) agreement for the V773 Tau A component mass range. The models 
indicate that the V773 Tau A components are still in their quasi-isothermal ("Hayashi") 
contraction phase. More quantitatively some model predictions seem to match the A com- 
ponent properties better than others. For instance, the solar-abundance Siess et al (2000) 
(upper left) and Palla & Stahler (1999) (lower-right) tracks near the dynamical masses tend 
to under-predict our estimated component temperatures and luminosities. In particular for 
the Siess et al (2000) models, this is similar to our previous results for HD 98800 B (Boden 
et al. 2005b). This apparent temperature discrepancy seems to reinforce a suggestion by 
Montalban et al (2004) that the solar-abundance Siess et al (2000) models appear to be sys- 
tematically cool compared with other models. Conversely, the D'Antona & Mazzitelli (1997) 
tracks (lower-left) at our mass values appear to slightly over-predict our component temper- 
atures and luminosity estimates. The Y 2 model predictions (upper-right; including tracks 
for the specific masses of the components) are seen to match our parameter determinations 
very well. 

Figure 6 gives a final model comparison. Recently Montalban & D'Antona (2006) have 
updated D'Antona & Mazzitelli (1997) models with a revised "2D radiative-hydrodynamic" 
(2D RHD) convection treatment. Figure 6 shows the result of this updated modeling with 
tracks from D'Antona & Mazzitelli (1997, in black) and Montalban & D'Antona (2006, 
in blue); the presentation is similar to that given in Figure 5, except we have displayed 
only mass tracks near the component dynamical mass estimates (1.3 and 1.5 M ). The 
2D RHD convective treatment increases the convective "efficiency" and energy transport in 
the atmosphere, moving the tracks toward lower temperatures during the "Hayashi" phase, 
and brings the model predictions into excellent agreement with the component parameters 
presented here. 

Isochrones from the models we considered would predict a range of component ages from 
2-7 Myrs, but restricting attention to the Montalban & D'Antona (2006) and Y 2 models, 
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Fig. 5. — V773 Tau A Components Compared With PMS Models. Here we show the 
V773 Tau A components in theoretical HR diagrams (luminosity/T e fj ), and solar-composition 
PMS evolutionary tracks and isochrones by Siess et al (2000) (top left), Y 2 (Demarque et al 
2004) (top right), D'Antona & Mazzitelli (1997) (bottom left), and Palla & Stahler (1999) 
(bottom right). All panels show available tracks in the 1.2 - 1.8 M range, and isochrones 
at 1, 3, 5, 7, and 10 Myrs. In the each panel we highlight the available tracks that best 
corresponding to our dynamical mass values; the Y 2 models uniquely allow us to compute 
tracks at these values. There is excellent agreement between our component parameter 
estimates and the predictions of the Y 2 tracks; other models included here are less consistent 
with the inferred component parameters. 
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a more plausible age estimate is 3 ± 1 Myrs, with both components appearing coeval to 
within present uncertainties. Such a system age is broadly consistent with age estimates 
for the Taurus population using a variety of evolutionary models (e.g. WG2001; Palla & 
Stahler 2002; Luhman 2004). More specifically this age estimate would put V773 Tau A in 
the putative intermediate-age population for Taurus as defined by Palla & Stahler (2002). 

4.2. Summary and Conclusions 

We have presented new interferometric (Tables 1 and 2) and RV data (Table 3) on the 
PMS binary subsystem V773 Tau A. Modeling these data in conjunction with older VLBI 
astrometry from P1996 allows us to estimate the visual and physical orbits for V773 Tau A, 
physical parameters of the components, and place significant constraints on PMS stellar 
evolution models. We characterize the V773 Tau A orbit model presented here as preliminary; 
we plan additional RV, VLBA, and KI observations that will further refine the physical orbit 
model for the A subsystem and constrain the component parameters. Further, ongoing RV 
and astrometric observations will continue to in order to probe the A-B orbit, and this in 
turn will refine the A subsystem RVs and orbit. Details of A-B orbit modeling will be 
documented in a future publication (Torres et al 2007). 

The V773 Tau physical parameters inferred from our orbit model are in good agreement 
with previous results. The most notable difference in our results with previous work is in 
the system distance; our distance estimate is at modest variance with the VLBA parallax 
from L1999. However, the L1999 data reduction did not account for the binarity of A 
which is clearly resolved in multiple epochs (P1996, this work). With that caveat we believe 
the L1999 distance and ours are in acceptable agreement. Because our distance estimate 
depends only on the apparent and physical orbits, we believe it is likely accurate at its 
stated uncertainty. We estimate our relative component mass errors at 9.1% and 7.3% for 
Aa and Ab respectively, making them comparable with other PMS stellar dynamical mass 
determinations (Hillenbrand & White 2004; Mathieu et al 2007). 

It is noteworthy that the P1996 and present VLBI data agree well with the A-orbit model 
presented here (including a model that excludes the VLBI data themselves, e.g. see Figure 2). 
Narrowly this would seem to support the P1996 interpretation that both A components are 
active radio emitters. More broadly, the ability to compare the radio morphology with 
at-epoch estimates of the component positions would seem to present an opportunity to 
constrain the radio emission and magnetic field topology in these stars. A more complete 
analysis of the VLBI data in this context will be made in a forthcoming publication (Massi 
et al 2007, in prep). 



- 20 - 



5 Klyr 3 Klyr l' Myr 
7 Myr 


11 


■ 10 Myr / 




1.5 M^ un / 1 
f if ^^^^ \ I i 


\ [ 


I if x / TI 
y I \ i-J — II- 

1.3 mjv \\, y j 


\ 1 \ 

1 X 
1 X 
1 \ 

IT v 


_ D'Antona & Mazzitelli 1997 Model* <\ \ \ 1 M > ' 
Montalban &DAntona 2006 Models " ; V v i\ \ 

V773 Tau Aa ■ — ' — ' \ \ x \ *\ 
V773 Tau Ah .— *— . 10Myr vVv W 

1 1 1 1 0, L x > 



7000 6500 6000 5500 5000 4500 4000 



Teff(K) 

Fig. 6. — V773 Tau A Components Compared With PMS Models From D'Antona & Mazz- 
itelli (1997) and Montalban & D'Antona (2006). The presentation is similar to Fig 5, but 
here we focus on tracks near the A-component dynamical masses (1.3 and 1.5 M ). The 
Montalban & D'Antona (2006, tracks in blue, isochrones in heavy dotted lines) models up- 
dated the older D'Antona & Mazzitelli (1997, tracks in black, isochrones in light dotted 
lines) models with a revised convection model, bringing the updated models into excellent 
agreement with the mass and radiometric parameters estimated here. 



- 21 - 



V773 Tau A radiometric estimates must necessarily account for (and are complicated 
by) extinction. Our estimate of the extinction is in reasonable agreement with earlier work 
(e.g. W1995; Ghez et al 1997; White & Ghez 2001), but prefers a slightly higher extinction 
when constrained by the spectroscopic temperatures (which are quasi-independent of the 
extinction). With this uncertainty in mind, we have made relatively conservative estimates 
of the errors in the radiometric properties. 

Even with the relatively conservative error estimates, the combination of our compo- 
nent dynamical mass and radiometric parameter estimates appear to differentiate between 
the range of solar-abundance PMS models considered here. In particular, we find better 
agreement between the estimated component parameters and predictions of solar-abundance 
models from Y 2 and Montalban & D'Antona (2006) than other models considered here, and 
these two models agree on an V773 Tau age estimate of approximately 3 ± 1 Myr; this is in 
good agreement with the multiple epochs of Taurus star formation (Palla & Stahler 2002). 
This age would place V773 Tau in the second of three age groups defined by Palla & Stahler 
(2002), i.e., stars with age 2 < t < 4 Myr. Such and younger stars are typically concentrated 
along the "central filaments" of the Taurus molecular clouds systems. V773 Tau appears in 
fact located at the periphery of a ~ 0.7 x 0.5 pc-sized "core" imaged in the high column 
density tracing C 18 molecule by Onishi et al. (1996, 1998, their core No. 5). From the 
observed intensity of the C 18 J = 1 — line one can derive an H 2 column density that 
translates into an Ay > 5. Extinction maps from Digitized Sky Survey source counts around 
V773 Tau's position also indicate an Ay around 5, but have a "resolution" of several arc 
minutes (Dobashi et al. 2005). These extinction values are several times larger than our 
and other extinction determinations for V773 Tau A (§ 3.1), and seem to indicate that the 
system is located just at the near border of a denser cloud core. 

Relative to the component parameters estimated here, the Siess et al (2000) and Palla 
& Stahler (1999) models predict cooler and less luminous components, while the D'Antona 
& Mazzitelli (1997) models predict hotter and more luminous components. However these 
conclusions would change if the component mass estimates were to systematically change 
(i.e. both components either more or less massive than present estimates) by only about two 
sigma. Qualitatively this seems relatively unlikely, but given the uncertainties in this study 
and the added complexity of multiple orbits in the V773 Tau system it would be prudent 
to withhold final judgment until the modeling of the A-B orbit improves with time and an 
expanded astrometric and radial velocity data set. 

Finally, it is important to note that placing constraints on stellar evolutionary mod- 
els in general, and PMS evolutionary models in particular requires precise measurement of 
more than just component dynamical mass, but also equally important parameters such as 
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luminosity, T c s, and abundance. Even at the present level of the work here (i.e. dynamical 
masses better than 10% 1-sigma), uncertainties in L and T e g- are significant factors in limiting 
conclusions we can draw about model predictions. As for abundance, we have restricted our 
comparisons of V773 Tau A component properties with solar-composition models. Funda- 
mentally this is a practical consideration; some of the model sets considered here only support 
solar abundances. Nominally abundances in Taurus- Auriga (and by consequence V773 Tau) 
are thought to be near solar values (see Padgett 1996). However, as in other sectors of 
the HR diagram, more secure abundance determinations in star forming regions, along with 
improved measurements of luminosity and temperature will be critical for definitive tests of 
PMS models. 
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